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Ground-state triply and doubly heavy baryons 
in a relativistic three-quark model 

A. P. MartjTienko* 

Samara State University, Pavlov Street 1, Samara 44^011, Russia 

Mass spectra of the ground-state baryons consisting of three or two heavy {b or 
c) and one Ught (u, d, s) quarks are calculated in the framework of the relativistic 
quark model and the hyperspherical expansion. The predictions of masses of the 
triply and doubly heavy baryons are obtained by employing the perturbation theory 
for the spin-independent and spin-dependent parts of the three-quark Hamiltonian. 
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The transition from two-quark bound states to three-quark bound states opens new prob- 
lems which refer both to the form of the quark interaction in the baryon and the structure 
of three-quark relativistic wave equation describing this system. In general form they were 
studied and solved already by many authors [l], UH, 0, H, 0, 0, 11] (see other references in 
Ref.jol). In practice it is important to have an approach which can allow to obtain simple 
and reliable estimates for the different experimental quantities regarding to the baryon spec- 
troscopy, the production and decay rates. Whereas heavy baryons with one or two heavy 
quarks were investigated both theoretically and experimentally, the triply heavy baryons 
fig-^QjQg containing b— and c— quarks have not studied so much. The estimate of masses 

Their 
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of the lowest-lying (ccc), (cc6), (bbc) and (bbb) states is presented in Refs.l 
production in a c— or b— quark fragmentation is calculated in Refs. [l^. [l3|. The doubly 
heavy baryons {Q1Q2I) represent a unique part of three-quark systems. Two heavy quarks 
compose a localized quark nucleus while the light quark moves around this color source 
at a distance of order (l/mg). This picture leads to the quark-diquark model for doubly 
heavy baryons which was used in Refs.[li,[ll,[i3 for the description of the mass spectrum 
and decay widths. Moreover, relativistic and bound state corrections to the mass spec- 
tra of mesons and baryons (in the quark-diquark approximation) and their different decay 
rates were also considered in the relativistic quark model 14, 15, 0, 17, 3]. Estimates for 
the masses of baryons containing two heavy quarks have been presented by many authors 



23| using different QCD inspired models for the quark interactions. The 



0, H H ii 

aim of the present paper is a twofold one. First, we go beyond the scope of the quark-diquark 
approximation in Refs. IJ, ll5| treating the total baryon Hamiltonian as a sum of two-quark 
interactions and using the hyperradial approximation for the ground state triply and dou- 
bly heavy baryons. Secondly, we take into account relativistic and bound state corrections 
of order v'^/c? by the perturbation theory. So, the purpose of our new investigation con- 
sists in the elaboration of an alternative calculational scheme of the baryon mass spectrum 
as compared with the earlier performed investigations in Refs. [3, 15] through the use of 
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the three-quark approach to the baryon problem formulated in Refs. [19|, |20|, |21[ with the 
Hamiltonian containing the spin- independent and spin-dependent corrections of order w^/c^. 

It should be notedjthat the theoretical results on the mass spectrum of triply and doubly 
heavy baryons [sl, 0, 0] remain untapped. In the past several years, the SELEX Collabo- 



ration has reported the first observation of doubly charmed baryons [2J, |25l]. But most 
recently, BaBar Collaboration has reported that they have not founded any evidence of dou- 
bly charmed baryons in e~^e~ annihilation j26|. Nevertheless, we can expect that the mass 
spectra and decay rates of triply and doubly heavy baryons will be measured before long. 
This gives additional grounds for new theoretical investigations of triply and doubly heavy 
baryon properties. 

In order to describe the mass spectra of baryons the different three-quark Hamiltonians 

H, El devoted to the 



are used 0, [3, H, 21 ■ The effective Hamiltonian of Refs.|l 
calculation of the mass spectra of doubly heavy baryons in the three-body approach is a 
sum of the string potential 1^™"^ and the Coulomb interaction potential V'-' . A consistent 
derivation of the three-quark potential was done in the Wilson-loop approach in Ref. [23, 2^ 
which accounts for corrections of order Accounting for the growth of the strong 

coupling constant for the interaction of a light quark with the heavy quark, we have included 
in the pure nonrelativistic three-quark Hamiltonian the vacuum polarization corrections of 



order in the form 



29, 30 



Ho 
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47r 



[ai + S^/ePo + 8/3o ln(/iij |ry|)) 



(1) 

We take the Hamiltonian (1) as the initial approximation in our study of the baryon mass 
spectrum. The operator of quark momenta Pi = —i^, fi is the position of quark i with 



respect to a common string-junction point which coincid es appr oximately with the center-of- 

Fj — Vj (compare 



mass point (for a more detailed discussion see Refs.HSlSl) and 

with FigJ_). The parameters of the confinement part of the potential are the following 



JJ,ll7|,ll8|: A 



0.18 Gel/^ B = 
0.33 GeV. 



—0.16 GeV . We take the quark masses rrih = 4.88 GeV , 
nis = 0.5 GeV as in our previous jH, 16, 3] calculations 



rric = 1.55 GeV, mu,d 

of different hadron properties on the basis of the relativistic quark model. For the dependence 
of the QCD coupling constant a^^ = cts{jl1j) on the renormalization point /2?,- we use the 



two-loop result 3ll 



as(/U, 



~2 ^ 



vr 



(3i InL 



/3o 



102 



38 



-n 



, L = ln(/i?./A 



(2) 

where A = 0.168 GeV, fiij = 2mimj/ {irti + irtj) is the renormalization scale, ai = (31 — 
10n//3)/3 and n/ is the number of flavours. The three-body Hamiltonian (1) and the 
Schrodinger equation can be reduced to the two-body form. For this aim, the three-body 



Jacobi coordinates become useful 



miTi + rrijV 



"J '-3 



nii + rrij 

where the coefficients aij, Pij are expressed in terms of the appropriate reduced masses: 




mij,k _ miTrij 



{rrii + mj)mk 
{rrii + rrij + ruk) ' 



(3) 
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/i is an arbitrary mass parameter which disappears in final expressions. Evidently, the 
coordinate p^j is proportional to the distance between quarks i and j, and the coordinate 
A^^- is proportional to the distance between the quark k and the center-of-mass of quarks 
Together with the center-of-mass coordinate Rem. the Jacobi coordinates determine 
completely the position of the system. 

In the center-of-mass frame (Rem. = 0) the operator of the kinetic energy can be written 
in the Jacobi coordinates p, A as follows: 

where K'^{Q) is the angular momentum operator, whose eigenfunctions (hyperspherical har- 
monics) are determined by the equation |34| : 

K\n)YKin) = -K{K + A)YKin). (6) 

Here fl designates five angular coordinates and R is the six-dimensional hyperradius 



R = sjpl + X% p = Rcos9, \ = Rsm9. (7) 
The baryon wave function \l'(p. A) can be presented as an expansion over functions Yx{^)'- 

^{p,x) = J2^K{R)YK{n). (8) 

K 




FIG. 1: The configuration of the three-quark system {QiQ2Q3)- O is the string-junction point. 
Qi,Q2 are the heavy quarks b or c. The quark is treated as a heavy quark 6, c in the triply 
heavy baryon or a Hght quark q in the doubly heavy baryon. 



Appearing here the radial wave functions \E'x(-R) satisfy to the system of differential 
equations jH, 3^]. However, in the further study of the ground state triply and doubly 
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heavy baryons we use the hyperradial approximation in which K = and the bound state 
wave function \1/ = \l/(-R) does not depend on the angular variables in the six- dimensional 
space. Using the Jacobi coordinates the Schrodinger equation for the three-quark system 
can now be transformed into the following form: 



5 d 



2fi[dR^^ RdRj 3§ Ip.^l 
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(9) 
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where coefficients •jij are given by 



ay a,- 



a,: 



fi{mi + rrij 



\ mfc(mi + m2 + m^) 



(10) 



Averaging Eq.(9) over angular variables, we can present the Schrodinger equation for the 
reduced radial wave function x(-R) = R^^^'^iR) in the form of a two-body wave equation: 



dR^ 



+ 2/X [E-V{R)]x{R)=0, 



V{R) = (^kR'' + ^2^, ao = ^B, ai 

k=-2 R ^ 



— AgT., a-2 = -, 
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The Schrodinger equation (11) determines the baryon mass spectrum in the initial ap- 
proximation. It can be solved numerically by the use the Mathematica program 35|]. The 
baryon wave function ^(2;) [x = R - ^JJi) for the state (ccc) is shown in Fig.2. 

The next step in the solution of the spectral problem is related to the consideration of 
different corrections to the Hamiltonian H^. In the present study we suggest that the three- 
quark potential has the form of the sum of twin quark the Breit-like interactions (we neglect 
orbital motion of the quarks) [H, [l^, 17, 18|: 



V^^^(ri, r2, rs) = V^^(ri, r2, rg) + V'^^x^, r^, rg) + ^ ^Mr,, r2, rg) 



(13) 
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FIG. 2: The wave function of the baryon (ccc) obtained after numerical solution of the Schrodinger 
equation (11). 
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where {...}vf denotes the Weyl ordering of operators. Accounting that the hyperradius R 
is independent of the order of the quark numbering we can average the potentials (14)-(19) 
over the functions ^(-R) to find their contributions to the baryon mass. The basic relation 
used for this aim is the following: 



^^{R)AVk{p,X)^{R). 



(20) 



Let us note, that the volume element dpdX can be written in terms of hyperradius R and 
the angle 9 {p ^ RcosO, A = i?sin^, < ^ < 7r/2): 



dpdX = {'in)^R^ sin^ cos^ 0dRd0. 



(21) 



Having the solution of Eq.(ll) in the numerical form, we have calculated all corrections 
from the interaction operators (14)-(19) numerically. In Table I we present the values of 
corresponding spin-independent matrix elements for the three-quark system {Q1Q2Q3) with 
the precision 0.001 GeV. This is not the accuracy of the mass spectrum calculation because 
the different theoretical uncertainties remain large (sec the discussion below). A number of 
the Breit-like potentials (14)-(19) is dependent on the relative quark distances r^j which could 
be expressed directly through the variables p^^ : r^j = p^Jotij. The confinement contributions 
and the pure relativistic corrections Z^f=iP|/8mf in the potential (13) are determined in 
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TABLE I: Contributions of the spin-independent terms of the potential (13) to the mass spectrum 
of the triply and doubly heavy baryons (in GeV). 



Quark 
content 


< AVi > 


< AV2 > 


< AV3 > 


< AV4 > 


< AV5 > 


< AVe > 


Summary 
contribution 




(ccc) 


0.032 


0.010 


-0.080 


-0.019 


-0.001 


0.001 


-0.057 




(ccb) 


0.016 


0.009 


-0.059 


-0.012 


-0.002 


0.0001 


-0.048 




(bbc) 


0.007 


0.007 


-0.033 


-0.008 


-0.002 


-0 0003 


-0.029 




(bbb) 


0.005 


0.003 


-0.018 


-0.006 


-0.001 


-0.001 


-0.018 




(ccq) 


0.037 


0.044 


-0.243 


-0.037 


0.020 


0.013 


-0.166 




(ccs) 


0.046 


0.025 


-0.167 


-0.029 


0.004 


0.006 


-0.115 




(bbq) 


-0.035 


0.039 


-0.131 


-0.012 


0.013 


0.003 


-0.123 




(bbs) 


-0.009 


0.022 


-0.083 


-0.010 


0.001 


0.001 


-0.078 




(bcq) 


-0.002 


0.042 


-0.214 


-0.023 


0.016 


0.008 


-0.173 




(bcs) 


0.015 


0.024 


-0.143 


-0.018 


0.002 


0.003 


-0.117 





terms of the variables A- : r^. = — 7jj A- . The perturbative and nonperturbative relativistic 



corrections are grouped together for the convenience in the expression AV3. The sign of the 
first term in the potential (14) differs from the corresponding expression in Ref. 27|]. The 
reason consists in the values of the universal Pauli interaction constant k = —1 and the 
mixing coefficient e = — 1 in the relativistic quark model fixed from the analysis of heavy 
quarkonium masses and radiative decays. 

Another part of the total Hamiltonian depends on the quark spins. Neglecting the quark 
orbital momentum we present the spin-dependent part of the potential in the simple form 



1^ 



\/^^(ri,r2,r3) = Ec^,S.S,5(|r,,|), 



drriimj 
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(22) 



(23) 



We assume that the color-magnetic interaction (22) can also be treated perturbatively as 
potentials (14)-(19). The matrix element of the operator (22) can be transformed as follows: 



SD 



"9" ^ ~ 



< SjSj > a 



rriimjTC 



X^dX\^{X)\^ X 



(24) 



X 



TT \ 3 3 \mi + rrij ^nii — rrtj 



1 m,- + 
+ ^— 




To obtain the contribution of the spin-dependent interaction (22) to the energy spectrum 
we have to perform the addition of the quark spins Sj to the baryon spin S. Designating the 
heavy quark spins Si and S2, we introduce the spin of a doubly heavy diquark S12 = Si + S2. 
Then, the baryon spin is S = S12 + S3. Considering that the baryon wave function '^Si2SS^ 
is the eigenfuction for the operators S^g; '^^5 ^z, we can express the scalar product (S1S2) 
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TABLE II: Mass spectrum of ground states of triply heavy baryons (in GeV). {Q1Q2} denotes 
the axial vector two-quark combination (diquark). The hyperfine splitting has been neglected in 
Refs.0, 



Baryon 


Quark content 


JP 


This work 


[9] 


flO] 






(ccc) 


3/2+ 


4.803 


4.79 


4.925 


4.76 


^ccb 


{cc}b 


1/2+ 


8.018 








0* 


{cc}b 


3/2+ 


8.025 


8.03 


8.200 


7.98 


^bbc 


{bb}c 


1/2+ 


11.280 










{bb}c 


3/2+ 


11.287 


11.20 


11.480 


11.19 


^bbb 


(bbb) 


3/2+ 


14.569 


14.30 


14.760 


14.37 



TABLE III: Mass spectrum of ground states of doubly heavy baryons (in GeV). {Q1Q2} denotes 
the axial vector two-quark combination (diquark), [Q1Q2] denotes the scalar one. 



Baryon 


Quark content 


JP 


This work 


m 


[22] 


[23] 


[4] 


[21] 




{cc}q 


1/2+ 


3.510 


3.620 


3.478 


3.66 


3.61 


3.69 


"cc 


{cc}q 


3/2+ 


3.548 


3.727 


3.61 


3.74 


3.68 






{cc}s 


1/2+ 


3.719 


3.778 


3.59 


3.74 


3.71 


3.86 


n* 


{cc}s 


3/2+ 


3.746 


3.872 


3.69 


3.82 


3.76 






{bb}q 


1/2+ 


10.130 


10.202 


10.093 


10.34 




10.16 


"bb 


{bb}q 


3/2+ 


10.144 


10.237 


10.133 


10.37 






^bh 


{bb}s 


1/2+ 


10.422 


10.359 


10.18 


10.37 




10.34 


^''bb 


{bb}s 


3/2+ 


10.432 


10.389 


10.20 


10.40 








{cb}q 


1/2+ 


6.792 


6.933 


6.82 


7.04 




6.96 


^'cb 


[cb]q 


1/2+ 


6.825 


6.963 


6.85 


6.99 






"cb 


{cb}q 


3/2+ 


6.827 


6.980 


6.90 


7.06 






^cb 


{cb}s 


1/2+ 


6.999 


7.088 


6.91 


7.09 




7.13 


Kb 


[cb]s 


1/2+ 


7.022 


7.116 


6.93 


7.06 






Kb 


{cb}s 


3/2+ 


7.024 


7.130 


6.99 


7.12 







through the diquark spin S12: 



S1S2 



5*12(5' 



12 



' 2. 



(25) 



where 5i2 = (the scalar diquark) or 5i2 = 1 (the vector diquark). For averaging the two 
other spin dependent terms (S1S3) and (S2S3) there is a need to use the unitary transfor- 
mation from the wave functions ^ Si2SSz to the eigenfunctions SnSS^ (or ^ S2zSS:,) i36j: 



>Sl3 



4 



i5i+52+S'3+S' 



v/(25i2 + l)(25 




53 5i 5 
5*9 5* 5* 



13 



SisSSz 



(26) 



12 
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Corresponding values of 6j - symbols are taken from Ref.[36j. The presence of operators 
(S1S3) and (S2S3) leads to a mixing between states with different spin 5*12 of the doubly 
heavy diquark and definite value S = 1/2. In the case of the (bcq) and (bcs) baryons we 
have the following mixing matrices: 

/ -0.005 0.015 X / -0.004 0.010 X 

V 0.015 - 0.020/^ ^' V 0-010 -0.015/^ ^ ^ ^ 

After the matrix diagonalization we obtain the numerical values for the masses of the 
baryons {bcq) and {bcs) = 1/2"'", 3/2"'" which are presented in Table III. In the case of 
triply and doubly heavy baryons with two identical heavy quarks c or b, the diquarks (cc) 
and {bb) have the spin 1 (the axial vector diquark). The hyperfine mass splittings for doubly 
heavy (cc) and (bb) baryons are the following: AM(Scc) = 0.038 GeV, AM{n^) = 0.027 
GeV, AM(Sbb) = 0.014 GeV, AM{Qbb) = 0.010 GeV. They are more than two times smaller 
as compared with the splittings obtained in Ref . fl3] . Whereas the masses of triply heavy 
baryons are in the agreement with the earlier performed calculations in Refs.jo], 10, 11| the 



obtained masses of doubly heavy baryons lie for the most part lower than our predictions 
from jr^. We expect that higher order corrections both in as and 1/m could change these 
results. 

Special attention has to be given to the accuracy of the performed calculation. In the 
case of triply heavy baryons the expansion over \pb,c\/^b.c is well defined. Numerical es- 
timate of relativistic effects gives the following expectation values: < pl/m^ >^ 0.34, 
< pl/ml >~ 0.09. So, the next to leading order contribution p^^/m^ leads to the theo- 
retical uncertainty ±0.030 GeV for the charm baryons and ±0.005 GeV for the bottom 
baryons (see Table II). The presence of the light quark in the baryon leads to the increase 
of the contribution of the relativistic effects connected with the motion of quarks q and s 
to the mass spectrum. Indeed, we can estimate the average value of the square momen- 
tum of light quark using the obtained baryon wave function (see Fig. 2). Introducing the 
parameter < p'^/m'l > where the brackets designate the matrix element with the func- 
tions of Eq.(9) we find that in the case of light quark q it has the value of order unity. 
So, summary perturbative and nonperturbative relativistic contribution is sufficiently large 
(see third column of Table I). Nevertheless, we suppose that three terms of the expansion 



^3 = \/P3±'^3 ~ ^3 ± P3/2^3 — Ps/Smg represent the light quark relativistic energy 
ea with sufficiently high accuracy (nearly 10%). The most significant error is connected 
with the relativistic corrections in the second order perturbation theory (PT). The reduced 
Green's function (^(R, R') for Eq.(ll) can be constructed numerically using the solutions 
of the Schrodinger equation (11) for different principal quantum numbers n. Saturating the 
sum over n in the expression of the Green's function by 10-15 excited states, we obtain ap- 
proximate value G(R,, R') which can be employed in the second order PT. Our preliminary 
estimate of the relativistic contributions in the second order PT shows that these corrections 
can comprise near 30% of the contribution < AV3 > (see Table I). 

Let us note that exists another possibility to rationalize the relativistic energy of the light 
quark. It demands the introduction of the effective energy (the momentum squared) of the 
light quark in the bound state. Indeed, we can use the following transformation of the light 
quark energy: €3 = ^3+^? = ± where = E3/2, and the light quark 

effective energy E3 can be expressed in terms of the baryon mass: E3 = Mb— mi —m2. In this 
case the quantity "fh^ plays the role of new effective mass of the light quark in the baryon and 
the addendum ml/ is equivalent to the term 777.3/2 used in Refs. 19, 2^, 21 1 if it is granted 
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that the effective mass tji^ coincides with the composite quark mass m^. Numerical value of 
the light quark effective energy E3 ^ 0.4-=-0.5 GeV, and the shift rr^ E2, ~ 0.21-=-0.22 GeV is 
in the agreement with the value of the term in Refs. [Tol. 20 . 21|. The appearance of the 
effective quantities in this approach evidently leads to the definite theoretical uncertainty in 
the baryon mass calculation. So, the presence of the light quark in the doubly heavy baryon 
essentially complicates the mass spectrum calculation with the sufficiently high accuracy. 
To improve the theoretical results of Table III the corrections of the second order PT should 
be calculated together with the matrix elements of the potential up to the order. 
The ground state wave function of the triply charm or bottom baryons is evidently the 
angle independent. The used hyperradial approximation is valid for it with high accuracy. 
For triply heavy baryons (cc6), {hhc) or for doubly heavy baryons {QiQ2(l) this hyperradial 
approximation is less applicable. So, the dependence of the baryon wave function on the 
angular variables in the six-dimensional space also has to be taken into account. 
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